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Abstract 

Perturbations around autonomous one-dimensional single- species reaction- 
diffusion systems are investigated. It is shown that the parameter space 
corresponding to the autonomous systems is divided into two parts: In 
one part, the system is stable against the perturbations, in the sense that 
largest relaxation time of the one-point functions changes continuously 
with perturbations. In the other part, however, the system is unstable 
against perturbations, so that any small perturbation drastically modifies 
the large-time behavior of the one-point functions. 
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1 Introduction 



Reaction-diffusion systems, is a well-studied area. People have studied reaction- 
diffusion systems, using analytical techniques, approximation methods, and sim- 
ulation. The approximation methods may be different in different dimensions, 
as for example the mean field techniques, good for high dimensions, generally 
do not give correct results for low-dimensional systems. A large fraction of an- 
alytical studies, belong to low-dimensional (specially one-dimensional) systems, 
as solving low-dimensional systems should in principle be easier. [1—13] . 

One of the reasons people want to find and solve exactly-solvable systems, is 
that one can use perturbative methods to investigate other systems, which are 
not exactly solvable but near some exactly-solvable systems. The term exactly- 
solvable have been used with different meanings. For example, in [14], [15], 
and [16], solvability (or integrability) means that the TV-particle conditional 
probabilities' S-matrix is factorized into a product of 2-particle S-matrices; while 
in [17-26], solvability means closedness of the evolution equation of the empty 
intervals (or their generalization). 

In [27], a ten-parameter family of reaction-diffusion processes was introduced 
for the systems among which, the evolution equation of n-point functions con- 
tains only n- or less- point functions. We call such systems autonomous. There, 
for these models the average particle-number in each site was obtained exactly. 
In [28,29], this has been generalized to multi-species systems and more-than- 
two-site interactions. 

Among the important aspects of reaction-diffusion systems, is the phase 
structure of the system. The static phase structure concerns with the time- 
independent profiles of the system, while the dynamical phase structure concerns 
with the evolution of the system, specially its relaxation behavior. In [30-33], 
the phase structures of some classes of single- or multiple-species autonomous 
reaction-diffusion systems have been investigated. These investigations were 
based on the one-point functions of the systems. 

In [34] , the two-point functions of autonomous single-species translationally- 
invariant one-dimensional reaction-diffusion systems were studied. The two- 
point function for such systems was obtained, and it was shown that it exhibits 
a non-trivial dynamical phase structure. The dynamical phase structure of the 
system was also investigated. 

In this article, we want to use perturbation to study systems which are 
not exactly autonomous, but nearly autonomous. By this, it is meant that 
the rates of these systems are equal to those of an autonomous system, plus 
a small perturbation. The scheme of the paper is the following. In section 
2, autonomous systems are briefly introduced. In section 3, non-autonomous 
perturbations around autonomous systems are considered, and their effect on 
the evolution of one-point functions is investigated. From this investigation, 
it turns out that some of the autonomous systems are unstable with respect 
to perturbations, in the sense that any small perturbation drastically modifies 
the large-time relaxation of the one-point functions. Section 4 is devoted to a 
concrete example. 
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2 Autonomous systems and the evolution equa- 
tions 



To fix notation, let's briefly introduce the autonomous systems. Consider a one- 
dimensional periodic lattice, every point of which either is empty or contains 
one particle. Let the lattice have L + l sites. The observables of such a system 
are the operators N", where i with 1 < i < L + 1 denotes the site number, 
and a = 0, 1 denotes the hole or the particle: is the hole (vacancy) number 
operator at site i, and is the particle number operator at site i. One has 
obviously the constraint 

s a N? = 1, (1) 

where s is a covector the components of which (s Q 's) are all equal to one. The 
constraint simply says that every site is either occupied by one particle or 
empty. A representation for these observables is 

N? := 1 ® ■ ■ ■ ® l ®N a ® 1 ® • • • ® 1 , (2) 

i—l L+l-i 

where N a is a diagonal 2x2 matrix the only nonzero element of which is the 
a'th diagonal element, and the operators 1 in the above expression are also 2x2 
matrices. It is seen that the constraint Q can be written as 

s-N=l, (3) 

where N is a vector the components of which are N°"s. The state of the system 
is characterized by a vector 

PeV<g> •••<&¥, (4) 



where V is a 2-dimensional vector space. All the elements of the vector P are 
nonnegative, and 

SP = 1. (5) 

Here S is the tensor-product of L + 1 covectors s. 

As the eigenvalues of the number operators N" are zero or one (and hence 
these operators are idempotent), the most general observable of such a system 
is the product of some of these number operators, or a sum of such terms. Also, 
the constraint shows that the two components of N, are not independent, 
so, one can express any function of in terms of 

ni := a • Nj, (6) 

where a is an arbitrary covector not parallel to s. Our aim is to study the 
evolution of the one-point functions ((nj's). 

The evolution of the state of the system is given by 

P = HP, (7) 
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where the Hamiltonian TL is stochastic, by which it is meant that its nondiagonal 
elements are nonnegative and 

SH = 0. (8) 

The interaction is nearest-neighbor, iff the Hamiltonian is of the form 

L + l 

H = J2 H i,i+u (9) 
t=i 

where 

Hj ti+ i := 1® ■ ■ •®l ®ff ® 1® - - -(a 1 . (10) 

i— 1 L — i 

(It has been assumed that the sites of the system are identical, that is, the 
system is translation- invariant. Otherwise H in the right-hand side of IjlUfl 
would depend on i.) The two-site Hamiltonian H is stochastic, that is, its non- 
diagonal elements are nonnegative, and the sum of the elements of each of its 
columns vanishes: 

(s®s)# = 0. (11) 

Using 

s ® s(a • N) <g> (b • N)H = a a bp H af3 jS s ® sTV 7 ® N s , (12) 

where a and b are arbitrary covectors, one can write down the evolution equa- 
tions of the one-, two-, or more- point functions of n^s. The evolution equation 
for the one-point function is 

^(n;> = a a sp H a ^ s (N? iVf +1 ) + s a a p H a ^ 5 (N^ Nf). (13) 

It is seen that the right-hand side of the above equation, contains two-point 
functions. In fact, in the evolution equation of n-point functions, there are 
generally up to (n + l)-point functions. However, there are systems for them in 
the evolution equation of n-point functions, only up to n-point functions arise. 
These are the autonomous systems. For a system with the Hamiltonian H° to 
be autonomous, following constraints hold among the reaction rates nd their 
reaction rates [27-29]. 

Q 7 sg + e 2 A a s s 7 , (14) 

where 

1 AOa rr°«/5 

A. 7 <5 --S/3 tl 7 s 

2 A 0a lS := S pH^ a lS . (15) 
It is not difficult to see that the constraints (|14fl are equivalent to 

ff°u(g)U = Au(giu, (16) 
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where 

u := (_\) , (17) 

and it is obvious that 

s-u = 0. (18) 

Now, consider an autonomous system satisfying the constraints l|14|l (or 
equivalently (|To)) ). and take the vector v satisfying 

E tA v =0, 

d,e-l J 

S-V=l, (19) 

and the covector a such that 

au=l, av = 0, (20) 

that is, the basis {a, s} is dual to {u, v}. In [28,29], it is shown that the matrix 
in the left-hand side of the first equation in i|19|) , has a left eigenvector with the 
eigenvalue zero. (This left eigenvector is s.) So it does have a right eigenvector 
with the eigenvalue zero as well. That is, there does exist a vector v satisfying 
<|19|) . In fact, one can even find a real vector v satisfying H19|) . From now on, a 
in © is assumed to satisfy <|2ll|) . 



3 Perturbations around autonomous systems 

Consider a system with the Hamiltonian H as 

H = H° + SH, (21) 

where H° is the Hamiltonian of an autonomous system, and the rates corre- 
sponding to 8H are small compared to those corresponding to H°. Our task is 
to investigate the evolution one-point functions corresponding to H , using per- 
turbation. As H is not necessarily autonomous, the evolution equation of the 
one-point function may contain two-point functions as well. However, the terms 
containing the two-point functions are proportional to the rates corresponding 
to SH , and hence are small. So, one can calculate the two-point function corre- 
sponding to the unperturbed system, and use it in the evolution equation of the 
one-point function of the perturbed system, to obtain up-to-first-order evolution 
of the one-point function of the perturbed system. 

3.1 The unperturbed solution 

Assuming that the initial condition is translationally-invariant, it is seen that 
the one-point function is independent of the site, and the two-point function 
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depends on only the difference of the sites' numbers. So, the evolution equation 
for the one-point function of the unperturbed system is 

^ = (/ + I /»)/°, (22) 



where 
and 



f ■= (ni) , (23) 



/i° =s <g> a H° u <g> v + a <g> s H° v <g> u, 

^° =s(g)aiJ v(g)U + a®siJ°u(g)v. (24) 

The one-point function /°, is easily seen to be 

f(t) = f(0)cM(v a + v°)t}. (25) 

Also, taking 

F?:=(n k n k+i )°, (26) 
(the two-point function of the unperturbed system) one arrives at 



dt 
dF? 
dt 



H°(F?_ 1+ F? +1 ) + 2v F?, 1 < i < L 
M F 2 + (^ + A K + p /° + a , (27) 



where 



A :=a®aif°u«u, 

=a <g> a i?° (u <g> v + v <g> u), 

=a<g>a#°v<g> v. (28) 
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It is seen that only five parameters enter the evolution equation of the up-to- 
two-point functions, and all of these can be expressed in terms of the matrix 
elements of 

3° ■= H° + ni?°n, (29) 

where II is the permutation matrix. These parameters can be rewritten as 





:=s <g> 




V 


v° 


:=s <g> 


ai?°v® 


u 


X° 


1 

:= 2 a 


® ai?°u 


® u 


P° 


:=a(g 


aH°u<% 


V 


a 


1 

:= 2 a 


® ai?°v 


® V. 



(30) 
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Taking a solution like 

i?(*)=£i? E o(0)exp(£ Q t), (31) 

E° 

it was shown in [34] that, the values of E° (energy-values) entering the two- 
point function are 0, := /j, + v , and any number in the interval Iq := 
[2v° - 2|m°|,22/° + 2|^°|], and possibly 

^ : = AO + ^ +A^o- (32) 

E% is among the possible values of E°, iff 

\v°\<\°-v . (33) 

The relation of E®, E%, and Jo, determines the relaxation behavior of the two- 
point function (its dynamical phase). Depending on the reaction rates several 
phases may occur [34]: 

I) Ei € Iq, and E® is not an energy. This is the slower phase, and the longest 
relaxation time is [— 2(v° — jit )] -1 . 

II) El € Iq, and E^ is an energy, in fact the largest one. This is the slowest 
phase, and the longest relaxation time is {— u° — A — [(/j,°) 2 /(A° — i/ )]} -1 . 

III) E® > Iq, and E\ is not an energy. This is the fastest phase, and the 
longest relaxation time is [— + /i )] -1 . 

IV) Ei > Iq, E§ is an energy, and E% < E®. This is the fast phase, and the 
longest relaxation time is [— + /z )] -1 . 

V) Ei > Iq, E® is an energy, and E% > E°. This is the slow phase, and the 
longest relaxation time is {— v° — A — [(fi°) 2 /(X° — z/ )]} -1 . 

3.2 The perturbed solution 

Now consider the Hamiltonian H, defined through J2J, which is not necessarily 
autonomous. Then defining e A a 7 a like Ijl5jl but with H instead of H°, it is 
seen that the evolution equation of the one-point function is 

d ^ 1 A a Jl\T7AT* \_i_„ 2 A a J AH AT$\ 



dt 



= a a 1 A a jS (N?Nt +1 ) + a a 2 A a ^(N^N?). (34) 



However, one cannot necessarily decompose e A a 7 s like (|14|) . Assuming trans- 
lational invariance of the initial conditions, one arrives at 

^ =a a { 1 A a 1 s + 2 A a lS )Ff, (35) 
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where 

/ :=(«-»), 

F? 5 :={N?Nf +l ). (36) 

As {u, v} is a basis, one can write F\ in terms of the tensor products of u and v. 
The corresponding coefficients can be found by multiplying the tensor products 
of a and s by F\. The result is 

= F lU ~<u 5 + f(u l v S + uV) + v 1 v 5 . (37) 

So, 

| = a a ( 1 A a 75 + 2 A a lS ) [F lU ^u s + /(«V + v~<u s ) + «V] . (38) 
Defining 

F> . — A 7( 5 A 7 <5, (*^9) 

l|3^|) recasts to 
df 

/l aOq. i 1 a Oct , 2 a Oct ( 2 ^Oq , 1 k?ck , 2 k?ck \ 

— =0 (l^ 7 s <5 + 2-4 5 S 7 + \A 7 Ss + 2 A JS 7 + £) ,j+ O ryg) 

x [f\uV + /(uV + uV) + vV] 
=(/i° + i/ )/ + ( 1 B 7<5 + 2 B l5 ) [F lU ->u s + /(uV + u V) + v^v s ] , (40) 

where 

e B l5 :=a a e B a lS . (41) 

As expected, the coefficients of F\ in the right-hand side are small (first order 
in terms of the perturbation SH). So, one can use the zeroth-order value of Fi 
in the right-hand side, to obtain the first-order value of /. 

From l|25[) and (|27[) . it is seen that if E® > Iq, then one can write 

*? = P ? + o /o° o /°» ( 42 ) 
/i u — A u — // u ^ 

where z satisfies 

/ i + // = / i°(z + z- 1 ) + 2^ , (43) 

and its modulus is less than 1, and F® is like (|31l) . but without a term corre- 
sponding to the energy- value E® = /jP + v° . So, one can write (|40|l like 

^ =(/x° + !/")/ + ( ^ 7i + 2 B l5 ) ( _£_ ^ z ulu& + u V + w V ) f 
+ { 1 B 1 s + 2 B lS ){F l u' 1 u s + «V). (44) 
This means that the energy- values entering /, are those entering i 7 ^ , and 

^ := ^ V-a' -m°/ *' (45) 
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where 



v 



=fjP + SfJL, 

=u° + Su, 



Sfi 

Su 



66 



=s (g &5H u (g v 
=s (g a Ji? v (g u 
=a (g s <5.ff u (g u 



(46) 



and 



SH := 5H + H5HU. 



(47) 



This shows that if E® is the largest nonzero energy- value entering F°, then E\ 
is the largest nonzero energy-value entering /. Otherwise, the largest nonzero 
energy- value entering / is the largest nonzero energy- value entering F®. So the 
relaxation behavior of /, can be deduced from that of F® as follows: 

I) In this phase the largest nonzero energy-value of / is 2(u° — fi°), and the 
perturbation causes a discontinuous change of the largest nonzero energy- 
value, from El = /i° + u° to 2(u° - 

II, V) In this phase the largest nonzero energy- value of / is E®, and the pertur- 
bation causes a discontinuous change of the largest nonzero energy- value, 



III, IV) In this phase the largest nonzero energy- value of / is E\, and the per- 
turbation causes a continuous change of the largest nonzero energy-value, 
from E® to E\. 

It is seen that the perturbation causes two different changes in the relaxation 
behavior of / (the one-point function). In the regions I, II, and V, the pertur- 
bation causes a discontinuous change in the relaxation behavior, which means 
that the autonomous system is unstable with respect to perturbations. In the 
regions III and IV, however, the relaxation behavior of / is continuous with 
respect to the perturbations, which means that the autonomous system is stable 
with respect to the perturbations, at least as long as first-order perturbations 
of the one-point function are considered. Mentioning one other thing is also 
in order: F® enters the evolution equation of /, iff 59 ^ 0. If 59 = 0, then / 
would contain only one (nonzero) energy-value, which is the one expected from 
changing one autonomous system to another. (Only ^ is replaced by \i and v° 
is replaced by v.) 

A real autonomous system, would in fact be only approximately autonomous. 
This means that there are always perturbations around the autonomous system. 
The above argument shows that only autonomous systems in the regions III 
and IV can be effectively autonomous. The parameter space corresponding to 
the effectively-autonomous systems is 



from El = fi° + i>° to E$. 




A° ^ (-//>°) - 1 + + 3(-m /^ )][1 - FgVg)T 



(48) 
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Figure 1 shows the regions corresponding to the autonomous system. (This is 
identical to figure 1 in [34].) 



4 An example 

Consider an Hamiltonian H° corresponding to an autonomous system: 



(- 


-3 + 


3lu 


OJ 


OJ 




1 - 


bJ 


— 2>oj 


oj 




1 - 


OJ 


OJ 


-3oj 




1 - 


OJ 


OJ 


OJ 



1 - OJ \ 

1-UJ 
1-LJ 

-3 + 3ujJ 



(-1 





\ 1 



The reactions of the corresponding system arc 



<DA 
A® 



any other state, 
any other state, 



with the rate w/4, 
with the rate w/4, 



AA 
AA 
AA 



%A, with the rate r + [(1 - w)/4], 

A®, with the rate r + [(1 - w)/4], 

AA, with the rate (1 - w)/4, 

<DA, with the rate (1 - w)/4, 

A0, with the rate (1 - w)/4, 

00, with the rate r + [(1 - w)/4]. 



For this Hamiltonian, one has 



and 



/x° = - 1 - 2r + 2w, 
v° = - 1 - 2r, 

p° =r. 



For the perturbation, consider the Hamiltonian 



SH = e 



(-1 o\ 





\l 0/ 



1 
1 

-v 



(49) 



(50) 



(51) 



(52) 



(53) 
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It is easy to see that this Hamiltonian does not correspond to an autonomous 
system. This Hamiltonian only increases the rate of the reaction AA — > 00 by 
e. Using this Hamiltonian, it is seen that 

S/j = —e, 
Su = —e, 

66 = -2e. (54) 

From l|52|) . it is seen that 

,,o 



= -1 



u° 1 + 2r 

A 1 



2 



(55) 



Comparing this with 148(1 . it is seen that the system is effectively autonomous, 
iff 

w> ^_^(l + 2 r). (56) 

As a special case of the above example, let us put u = 1. In this case, for 
the nonperturbed system we have the following reactions. 

%A — ► any other state, with the rate 1/4, 
A% — ► any other state, with the rate 1/4, 
00 -> 0A, with the rate r, 
00 -> A0, with the rate r, 

-> 00, with the rate r. (57) 

In this case, the system is effectively autonomous iff 

r<r := 1 - + ±. (58) 

It is seen that changing the value of r from to +oo, the system starts from 
phase III, passes through the phases IV, V, and II, and finally reaches the 
phase I. At r = r , the system goes from the phase IV to the phase V, which 
means that the system is no longer effectively autonomous. 
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Figure 1: the dynamical phase structure in the plane (— ^i°/^°, — X°/u°). The 
thick line is the boundary between effectively-autonomous systems (right side of 
the boundary), and effectively- nonautonomous systems (left side of the bound- 
ary). 



